ASYMPTOTICS OF EIGENFUNCTIONS ON PLANE DOMAINS 



DANIEL GRIESER AND DAVID JERISON 



Abstract. We consider a family of domains (£In)n>0 obtained by attaching 
an TV X 1 rectangle to a fixed set Qo = {(%, V) '■ < y < 1, —<t>{y) < x < 0}, for 
a Lipschitz function cj> > 0. We derive full asymptotic expansions, as TV — > oo, 
for the mth Dirichlet eigenvalue (for any fixed m £ N) and for the associated 
cigcnfunction on Qpf. The second term involves a scattering phase arising 
in the Dirichlet problem on the infinite domain f2oo. We determine the first 
variation of this scattering phase, with respect to rf>, at <j> = 0. This is then used 
to prove sharpness of results, obtained previously by the same authors, about 
the location of extrema and nodal line of eigenfunctions on convex domains. 
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1. Introduction 

For a Lipschitz function <j> ■ [0, 1] — ► [0, oo) and for TV 6 [0, oo] consider the plane 
domain (see Figure [lj 

(1.1) ^tv = {{x,y) e K 2 : < y < 1, -<t>{y) < x < N} 

and the eigenvalue problem for the Dirchlet Laplacian on fijv: 

(A + fi)u = on fijv 
u = at dil n- 

For TV < oo let //i(f2i\r) < ^(Qn) 5= • • • be the eigenvalues, counted with mul- 
tiplicities. The object of this paper is to study, for fixed m G N, the asymptotic 
behavior of ^m(^Jv) and of the associated eigenfunctions, as N — * oo. This will 
then be used to answer some questions left open in our study (see [8] HJH]) of the 
first and the second eigenfunction on general plane convex domains. 
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Figure 1 . The domain 
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Our first main theorem is: 

Theorem 1. There is a number a(</>) € [O,max0] such that for each m € N the 
with Dirichlet eigenvalue of Qn satisfies 



(1.2) 



yU m (f2 



N 



+ 0(N- 5 ), N 



oo. 



{N + a^)) 2 

In particular, the eigenvalues fii, . . . , u m of £1n are simple for N sufficiently large. 
The suitably rescaled eigenfunction u m satisfies, for all multiindices a, 



(1.3) sup \D a 

x>3 log N 

(1.4) 



t (x,y) - sin 



a(4>) 



sup 

ic<3 log N 



N + a(<t>) 
u m {x,y)\=0(N- 1 logN) 



sin Try | = 0(/V" J ), 



The constants in the error terms only depend on k, a and max <p. 

Thus, the spectral data are very close to the data obtained on the rectangle 
[— a(4>),N] x [0,1]. In fact, we will get complete asymptotic expansions for the 
eigenvalue and the eigenfunction and also much more precise information about the 
eigenfunction for small values of x, see Remark 1 121 and Theorem [15] (in connection 

with dsn])). 

The number a((j>) is closely related to a scattering phase in an associated non- 
compact problem, see Section [5] in particular Remark [5] Therefore, its dependence 
on <j) is very subtle. Our second main theorem gives a perturbation analysis of a(cf>) 
around = 0: 

Theorem 2. Fix a Lipschitz function <f) : [0, 1] — > [0, oo). Then, as e — ► 0, 



(1.5) 



a(e(f ) ) = 2e 0{y)(simry) 2 dy + 0(e 2 ). 
Jo 



As already mentioned, one motivation for the present study is to complement the 
results about the first and second eigenfunction u\, U2 on a plane convex domain 
Q, obtained in [8] [5]. In these papers we considered the maximum of u\ (which 
is w.l.o.g. assumed positive) and the nodal line of 112, and a central goal was to 
localize these objects in terms of corresponding objects for eigenfunctions of an 
associated ordinary differential operator. To state this more precisely, we first 
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normalize by a rotation and a dilation so that among all projections of Q onto 
lines the projection onto the y-axis has shortest possible length and this length 
equals one. Let n(x, y) = x be the projection map to the x axis. Let I = 7r(f2) and 
h : I — > (0,oo) be the 'height function' of il, that is, 

h(x) = length of the interval 7r -1 ({x}) n J7. 

d 2 

Let 4>\ , 4>2 be the first and second eigenfunction of the Schrodinger operator — — - + 

ax 2 

on /, with Dirichlet boundary conditions (defined in terms of the variational 

h(xY 

principle on Hq(I)). 

Theorem 3 ([8, 5]). Let the domain O be normalized as above and let h : I — > 
(0, 00) be its height function. Denote by M. the set where u\ achieves its maximum, 
by M — u^" 1 (0) the nodal line, and let {m} and {n} be the corresponding sets for 

01,02- 

There is an absolute constant C so that 

(1.6) ir(M) C [m-C.m + C], n(N) C [n - C, n + C]. 

Actually, Ai consists of a single point, by a well-known argument using the 
convexity of f2. Also, the uniqueness of the point m follows by standard arguments 
from the convexity of h, while the uniqueness of n is a general fact from Sturm 
theory. 

A consequence of the theorem is that, while clearly the interval / and the function 
h do not determine uniquely, these data do determine the location of the distinc- 
tive features of u up to a bounded error, uniformly for all domains normalized 
as above (in particular, uniformly as \I\ — > 00). 

The question left open in [8,5 is whether the result (|1.6p is sharp in order of 
magnitude as |/| — > 00. We will derive from Theorems [1] and [2] that this is in fact 
true: 

Theorem 4. There is c > and, for each N' > 0, a pair of domains fl, fl, 
normalized as above and with tt(Q) = 7r(0) of length N' and with the same height 
function, such that 

(1.7) dist(7r(M), n(M)) > c, dist(7r(A0, tt(AO) > c. 

See Figure [2] at the end of the paper. This should be put in contrast with the 
main result of [4 which states that the length of the interval tt(AT) is not bounded 
away from zero, but actually bounded above by C/\I\, for an absolute constant C. 

Our approach to Theorem[T]is via the method of matched asymptotic expansions. 
This is carried out in Section [3J Here one has to deal with two limiting problems: 
The first one is the equation {—j^i + \)ip(£) — on the interval [0, 1], and is easily 
solved explicitly. The second one is the equation (A + ir 2 )U = on the unbounded 
domain f^; we need some results about generalized eigenfunctions on fioo, which 
are known from scattering theory. For completeness we give a direct derivation of 
what we need in Section [21 Here, the quantity a(0) arises. Theorem [5] is proved in 
Section S] and Theorem [3] in Section 



4 



DANIEL GRIESER. AND DAVID JERISON 



The asymptotic behavior of spectral quantities on degenerating spaces similar to 
the family (Qn)n>o has been studied by many authors in different contexts. Re- 
garding invariants involving all the eigenvalues (like the determinant of the Lapla- 
cian) we only mention [7], [llj and [12] . Our results are also related to the inves- 
tigations of so-called thick graphs (or graph-like thin manifolds): When rescaling 
by a factor 1/N, one obtains a domain which is a l/2iV-neighborhood of a 
unit interval, except for a fixed (scaled) perturbation at the ends of the interval. 
Instead of the interval (considered as a graph with two nodes and one edge of length 
one connecting them) one may consider more general embedded graphs, and their 
e-neighborhoods (or more general 'e-thin' manifolds modeled on the graph). The 
convergence of eigenvalues, as e tends to zero, to spectral data on the graph it- 
self (then sometimes called a 'quantum graph') was studied in EH [T] , for the 
Neumann and closed problems. The Dirichlet and mixed boundary value problems 
are studied in the preprints [3] (where it is also proved that the asymptotic series 
constructed in this paper converge) and [10] , by different methods than the one 
used in this paper. The Dirichlet problem is more difficult to handle since the de- 
pendence on the counting parameter k appears in a lower order term, cf. equation 
(jl.2|) . While in these papers the graph edges are always straight lines, the case of 
a curved line (but without the perturbation at the end) is considered in [2J, where 
a nodal line theorem is proved. 



In this section we prove the results about generalized eigenfunctions on which 
will be needed in the proof of Theorem [TJ Since, for any m, the mth eigenvalue 
on f2jv converges to it 2 as N — > oo (by domain comparison), we need to consider 
the spectral value 7r 2 , which is the bottom of the continuous spectrum of — Aoir on 

Proposition 5. 

(1) There is a unique function U on fioo satisfying 

(A + 7r 2 )[/ = on 



2. Eigenfunctions on the infinite domain 



(2.1) 



U = at da 



1 oc 



U — x sin ny is bounded. 



(2) For this function U define 



(2.2) 




Then 



(2.3) 



U(x,y) 



{x + a((j))) sin iry + r{x, y) 



where the remainder decays exponentially as x — > oo, more precisely, 



oo 



(2.4) 



r{x,y) = ^r W (x)sinfc7ry, \d p x r {k) {x)\ < C v e~ kxl2 Vp,k. 



k=2 

Here, C p is independent of k and x and is bounded in terms o/max(/>. 
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Note that in the special case where cf> is a constant A, the function U is simply (x+ 
A)smiry, so a(4>) — A. Therefore, for general <p, the number a((f) tells how much 
the 'standard' problem, with fl^ = [0, oo) x [0, 1], has to be shifted so that its first 
(generalized) eigenfunction coincides asymptotically with that of the 'perturbed' 
problem. See also Remark [8] 

For simplicity, we assume all functions are real valued. Basic to all considerations 
is the explicit form of solutions of the homogeneous equation: 

Lemma 6. Assume u solves (A + ir 2 )u = in x > Xq, y € (0, 1), vanishes for 
y = and y = 1 and has at most polynomial growth as x — > oo. Then 

oo 

(2.5) u{x,y) = (Ai +Bix) shiny + A k e~^ k ' 2 ~ l7TX sinkiry forx>x , 



k=2 



for certain numbers A^, B\ 



Proof, u is smooth, so for each fixed x it is the sum of its Fourier series, u(x, y) = 
Sfcli u k(x) sin kiry, where Uk(x) = 2 J Q u(x, y) sin kiry dy. From (A + ir 2 )u = 
one gets u'£ + (1 — k 2 )-K 2 Uk = 0. This gives u\(x) = A\ + B\x and u^(x) = 
Ake~^ k2 ~ llTX + Bke^ k2 ~ l7rx for k > 2. Since u is polynomially bounded, so is Uk, 
and therefore Bk = for fc > 2. □ 

Lemma 7. Let w 6 L 2 (rioo) be supported in Q,\. 
(1) If u e -ff^fioo) so/ues 

(A + ir 2 )u = w on £1 
(2.6) u = at dn 



OO } 
OO ) 



U IS 



bounded 



then 



(2-7) hf H i {ni) < -C(u,w) 

where C is bounded in terms of the maximum of (p. 
(2) Problem (|2.6p has a unique solution u. 

Proof. (1) We integrate by parts and use the support assumption to obtain 

|2 J2 / „,2 



(2.8) -(u,w) = - u(A + Tr 2 )u = - uu n + \\7u\ 2 - n 2 u 2 

where u n is the outward normal derivative. Therefore, it is sufficient to prove the 
following two facts: 

r Ou(xq^ y) 

(2.9) If u is as in (12.51) and bounded then / u(xo,y) — - — dy < 0. 

Jo ox 

(2.10) 

IfuGF^fii) andu = 0onani\{x = l} then / \Vu\ 2 -tt 2 ( u 2 > c\\u\\ 2 h1 , Qi) 

./<> ./'.> 

To prove (|2.9[) . observe that B\ = in (|2.5p since u is bounded. Therefore, 
f Q u(x ,y)d x u(xo,y) dy = -7r/2 2^L 2 ^it^ 2 ~ lexp(-2\/fc 2 - lvrxo) < 0. 

To prove (|2.10p , consider the domain f2 which is the union of Cl±, the boundary 
piece {1} x (0, 1), and the reflection of Qi across this boundary. Since f2 C (—A, A) x 
(0, 1) for some A, the first Dirichlet eigenvalue of f2 is strictly bigger than tt 2 , so it 
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equals (it 2 + c)/(l — c) for some c > 0. The function on which equals u for x < 1 
and is symmetric with respect to the line {x = 1} is in _ffp(f2), so we can use it as 
test function and obtain (1 — c) f n |Vm| 2 > (it 2 + c) J n u 2 , which implies (|2.10[) . 
In this proof A and therefore c only depend on max <f>. 

(2) Uniqueness is clear from 1. To prove existence, we reduce to a compact 
problem. Define 'Dirichlet-to-Neumann operators' A^, A^, acting on functions on 
S = {1} x (0, 1), as follows: Given / € Hi (S) nH 2 (S), let vi,Voo be the solutions of 
(A + ir 2 )v = onf2i, D = (1 , oo) x (0, 1) respectively, with boundary values / at S 
and zero elsewhere, and i>oo bounded. Existence and uniqueness follows for v\ from 
the fact that the first Dirichlet eigenvalue of fli is bigger than ir 2 , and for Voo by 
explicit computation as in LemmaEl Set Nif := dvi/dx\s and -/V^/ := dvoo/dx\s- 
The restrictions exist and are in H 1 (S) since v\, Woo are in H 5 ^ 2 near S by standard 
regularity theory. 

From applied to u = we have (N^f, f) < 0, and from (|2~TU|) (Nif, f) = 
J dni M nVl =ifi 1 l Vwi l 2 - n v l ^ c II w iIIhMOi) - H^llffV2 C s). so we obtain 

(2.11) ((ATr - AU/J) > cll/ll^s), feHl(S)nH 2 (S). 

Along the same lines one sees that (Nif,g) < \\f\\H 1 / 2 (s)\\9\\H 1 / 2 (s) f° r f>9 G 
^□^(S 1 ), and this shows that A^i can be extended to a bounded operator Hq 2 (S) — > 
_ff _1 / 2 (S l ), and similarly for A^. By approximation, (|2.11| continues to hold for 

1 11 

/ 6 Hq' (S). This shows that Ni — Noo has closed range and therefore is surjective, 

1/2 ________ 

for if g G iJ ' (S 1 ) is orthogonal to the range then applying (|2.1ip to g implies g = 0. 

Now we find a solution u of (|2.6[) as follows: Let u be the unique solution of 
(A + tt 2 )v = w on Qi, viqq 1 = 0, and define / on S by (Ni — A^)/ = —dv/dx\s- 
For / find functions Vi, Voo as above. Define the function u by v + vi on and by 
Woo on D; at S, u has the value / from the left and right, and du/dx is dv/dx + Nif 
from the left and A^/ from the right, which are equal by construction. Therefore, 
u is the desired solution. □ 

Proof of Proposition^ (1) Let ip( x ) be a smooth function vanishing for x < 
and equal to one for a; > 1. Let u be the solution of (|2.6p with u> = — (A + 
tt 2 )(4>{x)x sin 7ry). Then [/ = it + i^(x)x sin 7ry solves (12.11) . From (|2 . T[) we have 
IMIff^ni) < C|l w IU 2 (ni) < C and therefore also 

(2-12) ||C/|| HI(ni) < C, 

with C only depending on max0. Uniqueness is clear from the uniqueness for u. 

(2) Apply Lemma [6] to U with xq = 0. Since U — x shiny is bounded, we have 
Bi = 1. Since, from (|2.12| and the trace theorem, J2T=i^t = \ So t^(0, t/) 2 c?j/ < 
C W U \\ 2 H 1 (n 1 ) - C '' tms S ives and (f2l| with a(0) = Ai, and evaluating 

Jo 1 [7(0, y) simrydy yields (g^l). □ 

Remark 8. We explain the relation of a(</>) to the scattering phase. Standard 
scattering theory (see, for example, [S]) yields that for s > close to zero the 
equation (A + tt 2 + s 2 )E = 0, i^an^ = has a unique polynomially bounded 
solution on of the form 



E s = [e- lsx + S( S y sx ) sin Try + r s {x, y) 
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for some number S(s) (the scattering matrix) and a remainder r s (x,y) of the form 
(|2.4p . The function S extends holomorphically to a neighborhood of zero in C and 
is real and of modulus one for real argument, hence may be written S(s) = e 17 ^ 
for a holomorphic function 7, the scattering phase. r s is also holomorphic in s and 
the estimates (|2.4p are uniform in s near zero, so one can take the limit s — > 
to get a solution of (A + tt 2 )U — 0. This solution is bounded, hence constant 
equal to zero by Proposition [5] Therefore S(0) = —1 and ro = 0. One can get 
a nontrivial solution by taking lim^o -E s = 4-, _ n E s , and this has leading term 

S CIS \s — 1. J 

— 2i(x + y(0)/2) sin7ry. Comparison with (|2.3[) then yields 
(2.13) a(4>) = ^(0). 

We will also need an extension of Proposition [5] 

Lemma 9. Assume v is a smooth function on fioo which vanishes at dQoo and for 
x > has the form 

v(x,y) =p(x) sin Try + r(x,y), 
with p a polynomial and r satisfying the estimates (|2.4ll . 
Then any polynomially bounded solution of the problem 

(A + tt 2 )u = v on f2oo 
(2-14) K ' 

U = at ofloc 

has the same form u(x,y) — q(x) sin7ry + s{x,y), where s satisfies the estimates 
HUD and 

(2.15) q"=p. 

Such solutions exist, and are unique up to adding multiples ofU, where U is defined 
in Proposition^ 

Proof. Let u be a solution of (|2.14[) . Taking the Fourier decomposition of u(x, ■) we 
get u(x,y) — q(x) sin7ry + ^2^=2 s ^ k H x ) sink-try, and then (|2.14p gives (|2.15|) and, 
for each k > 2, 

(2.16) (4 + 1 - fc 2 )s (fe) =r (k \ 

For any initial condition s^(0) — (|2 . 1 6[) has the unique polynomially bounded 
solution 

1 r°o 

(2.17) s {k) (x) =a k e~ ax / ( e ~ alx ~ zl - e - a{x+z) ) r {k) (z) dz 

where a — \/k 2 — 1. Since u is given and s^ k \x) — 2 J Q u(x, y) smknydy, we have 
a k = 2 Jq 1 u(0,y) sin kny dy. An easy calculation shows \dP.s^(x)\ < Ce~ kx l 2 , and 
this proves the first claim. 

To prove existence of a solution, choose qo satisfying q' ' — p, and a cutoff 
function x, equal to zero for x < 1 and to one for x > 2, and set uo(x,y) = 
x(x) (^qo(x) sin ny + ^2^ =2 s k \x) sin kiry^j. This is polynomially bounded. Then 

u = uq + h solves (I2.14p iff h solves (A + 7r 2 )h — v\, h^^ — 0, where v\ = 
v — (A + 7r 2 )uo is compactly supported. Lemma [7] guarantees the existence of a 
bounded solution h. 

Uniqueness is clear from Lemma [5] and Proposition [5] □ 
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3. Asymptotic expansions of eigenfunctions and eigenvalues 

In order to prove Theorem Q] we use the idea of matched asymptotic expansions. 
The strategy is this: First, we make reasonable guesses about the asymptotic be- 
havior, as N — > oo, of the mth eigenvalue and of certain scaled limits of the eigen- 
function. This leads to an ansatz in the form of formal asyptotic series in terms of 
powers of 1/N, whose coefficients are undetermined numbers (for the eigenvalue) 
resp. functions (for the eigenfunction) . The eigenvalue equation, the boundary con- 
dition and the condition that the various scaled limits must fit together ('match') in 
the transition region between different scaling regimes, yield a recursive system of 
equations for these coefficients. This system has a unique solution ( Proposition ! 1 Op . 
Given any approximation order, one then obtains a candidate for an approximate 
eigenvalue (by truncating the formal series), and also for an approximate eigen- 
function, which is obtained by a suitable patching of the data from the different 
scaling regimes. These candidates satisfy the eigenvalue equation with a small error 
(Proposition [13]), and from this we derive that they are close to actual eigenvalues. 
A domain comparison yields an a priori estimate on the actual eigenvalues, and this 
allows to conclude that all actual eigenvalues are obtained in this way, as well as a 
lower bound on the spectral gap. The spectral gap estimate then implies that the 
approximate eigenfunctions are close to actual eigenfunctions (Theorem ll5[) . Using 
this explicit information, it is easy to derive Theorem [T] 

3.1. The ansatz, formal eigenvalue and eigenfunction. Fix an integer m > 1. 
We want to find the mth eigenvalue and eigenfunction of f2jv, asymptotically as 
N — ► oo. In this and the next subsection we simply write A for the mth eigenvalue 
and u for an associated eigenfunction. As a guide, recall that for the unperturbed 
case fl N = [0, N] x [0, 1], we have (for N > m/2) 

A = 7r 2 + N m tt , u{x, y) — sin(N~ 1 mnx) sin(7ry). 

Our ansatz is guided by the following expectations: 

(1) The eigenvalue should have complete asymptotics: 

3A, 4 E R, i = 2,3,4,..., 

OO 

(3.1) A~7r 2 + ^iV _i Ai, N^oo. 

i=2 

Note that A = it 2 + 0(N~ 2 ) follows from domain comparison. 

(2) At any fixed x,y, suitably normalized eigenfunctions should converge as 
N — > oo, and even have complete asymptotics: 

Bfi-.Q^^R, i = 0,1,2,..., 

OO 

(3.2) u(x,y)~J2 N ~ l fi( x >y^ N^oo, (x,y) e fl^. 

i=0 

(3) When fixing £ = x/N and y 7 and letting N — > oo, u should converge, and 
even have complete asymptotics: 

3 9j : [0,1] x [0,1] -fR, j = 0,1,2,..., 

OO 

(3.3) u(m,y)~J2 N ~ j 9i&y)> N ^ °°> «,») G [0,1] x [0,1]- 
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We get conditions on all the coefficients from three sources: 

(I) The equation (A + X)u = 0. Formally inserting the asymptotics above, 
differentiating term by term, and successively equating powers of N, we 
get, for i, j = 0,1,2,...: 

i 

(f) (A+TT 2 )f i = -J2*kfi-k 

k=2 

j 

(g) (dy + TT 2 )g 3 = -<9f .9j-2 - x k9j-k 

k=2 

(terms with negative indices are set equal to zero). 
(II) Boundary conditions on dil n- (|3.2p . (|3.3|) give 

(bd f) fi = at arjoo 

(bdg) & = 0at{y = 0}u{y = l}U{£=l}. 

We will prove below that (f), (bd f) imply that each has the form 

i i 

(3.4) fi{x,y) = (p i (x)smivy + r i (x,y), <pj(x) = y^ J w a — ) 

1=0 

with ri satisfying condition (|2.4p . 
(Ill) Matching conditions. To ensure small errors when patching the fi and the 
gj to get an approximate eigenfunction, we need to correlate the large x 
behavior of the fi with the behavior at £ — of the gj . This is done by for- 
mally writing £\ N~ % fi(x, y) = £V N~ j gj(x/N, y), expanding ft according 
to (|3.4p and gj in Taylor series at £ = and equating the coefficients of 
N~*x . This gives Wu sin7r?/ = 9|^_;(0, y). This suggests to seek gj in the 
form 

(3.5) g j (€,y)=ip j (t)smny, 
and then the matching conditions read 

{M a ) ^ = ^d(o), l< t . 

Let us call a pair (J2i=o N~' i fi,Yl'j=o ^ 01 f° rma l series, with fi, gj of the 
form (|3.4p . (|3.5p . a formal eigenfunction with formal eigenvalue it 2 + Y2°Z. 2 if 
(f), (g), (bd f), (bd g) and (M^) are satisfied for all indices, and not both f , g are 
identically zero. Clearly, multiplying a formal eigenfunction by a non-zero scalar, 
that is a series N~ l a,i with ao ^ 0, yields a formal eigenfunction again, with 

the same formal eigenvalue. 

Proposition 10. If n 2 + J2i=2 N * s a formal eigenvalue then 

(3.6) A2 = m 2 7r 2 , for some meN. 

Conversely, for each m £ N i/iere is a unique formal eigenvalue with A2 = m 2 ir 2 , 
and the formal eigenfunction is unique up to multiplication by scalars. 
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Furthermore, we have 

oo 

\ " N-i \ ; = 

(N + af 



(3-7) 1^ N - 3 ^ = 7WT^ + 0{N -% 



j=2 

(3.8) ^Jv-^-(-)= S mm7r|^ + 0(iV- 3 ), 



where a = a(</>) is defined in 

The Xj, the coefficients of the (fj, and the constants in the estimates (|2.4|) of 
the remainders rj and in (|3 - T[) and (|3.8[) are all bounded in terms of j and the 
maximum of <f>. 

Proof. Because of Q3.5P we may rewrite (g), (bd g) as 

i+2 

WO (4 + A a )V»j = - X! A/V^+2-/ , 

1=3 

(bdV) ^-(i) = o, 

where we shifted the index by two. 

To prove (|3.6[) . note that fo is bounded and satisfies (A + 7r 2 )/o = on fioo, 
hence is zero by Lemma[7]j2). Then, (Moo) gives tpo{0) = 0, so we have 

(4 + A 2 )^ = 0, MO) - Ml) = 0, V^o # 0, 

and this implies (|3.6p . 

Now fix m € N. We construct a formal eigenvalue and formal eigenfunction with 
A2 = m 2 7r 2 , and satisfying the normalization condition 

(3.9) 2/ ^smm^d^Soj, j = 0,l,2,..., 

Jo 

and simultaneously prove its uniqueness. Since multiplying any formal eigenfunc- 
tion by the scalar ^2 Y^jLa J ^~ J Jo ^AO sinm7r£ d^j yields a formal eigenfunc- 
tion satisfying (|3.9p . this will prove the Proposition. 

First, by the argument proving ()3.6|) . and by (13.9j) . we must have 

fo = 0, ipo(0 = sinrnTT^. 
We now apply iteratively the following lemma. 

Lemma 11. Let J > 1. Given ip , . . .,ipj-i, fo, . . . A 2 , . . . ,Xj+i satisfying 

the equations (ip), (f), (Mji) and the boundary conditions for j < J, there are 
unique ipj, fj, \j+2 satisfying these equations for j — J and the normalization 

(EH). 



Proof. First, we choose a solution fj of (f) (with i = J) of the form (|3.4p . accord- 
ing to Lemma |9l removing the indeterminacy by prescribing wji — ?///_! (0). This 
determines wjq, and therefore ipj(0) by Mjq. Next, equation (■0) with j — J has a 
solution with given values at and 1, if and only if the right hand side satisfies one 
linear condition, and then the solution is unique up to multiples of ipo = sinrmry, 
therefore uniquely determined by condition (13. 9j) . The solvability condition is ob- 
tained by taking the scalar product of both sides with ipo and integrating by parts 
on the left. This gives 

(3.10) Vj(0HM0) = -A/ +2 /2, 
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where we have used (|3.9[) . and this determines Xj+2- 

To finish the proof, we only need to check that (Mji) is satisfied for I > 2. 
Now from (|2.15p . the polynomial ipj occuring in fj satisfies tp'j = — J2k=2 ^kfj-k- 
Equating coefficients of j^yX 1 ^ 2 simry we get wji = — ^2 k=2 Using 
the matching conditions (with j < J — 2) on the right and then the (I — 2)th 
derivative of equation (ip) , with j = J— I, we see that this sum equals d l ^ipj-i(0). □ 

The boundedness of all quantities in terms of max0 is also proved inductively, 
using the corresponding claims in Proposition [5] and Lemma [7J 

As an illustration, we carry this out for J — I: ipo{t;) = sinm7r£ gives W\\ = 
V'o(O) = m7r, this determines f\ = mnU, and (|2.3p yields wio = mna, hence 
ipi(0) = mira. Equation (|3.10|) now gives A3 = — 2m 2 7r 2 a, and we have — 
mira(l — ^) cos rrnra. 

It remains to check ()3.7|) and ()3.8() . The calculations of higher order terms can 
be simplified by introducing new variables x = x + a, N = N + a, expressing all 
functions in terms of £ and £ = x/N, and using formal series in N. This must give 
the same result (after changing variables and up to normalization) by uniqueness. 
We get in the J = 1 step: ipo{0 — sinm7r£ implies iin = mir, so 

fi = mirU 

(where U{x,y) = U(x,y)) as before, but now U{x) = x smiry + 0(e~ x ), so wio = 
and thus ^i(O) = 0, from which we get, using (|3.10p 

^1 =0,A 3 = 0. 

The J = 2 step yields W21 = 0, and since (A + t^ 2 )H = 0, we get from Proposition 
[5]that f 2 = 0. As before, this gives ip 2 = 0, A 4 = 0, and this proves (|3.7p . (|3.8p . □ 

Remark 12. The next term can be obtained as follows. For J = 3 we get W31 = 0, 
but now (A + 7r 2 )/3 = — A2/1 = —m 3 ir 3 U. To obtain an expression for 11130, we 
multiply this with U, integrate over Qa, apply Green's formula and let A — > 00. A 
short calculation (using (|2.15p also) gives 

w 30 = 2to 3 7t 3 6, b := lim ( [ U 2 dxdy - — ^ , 
A -^°° \Jn A 6 / 

and then (|3 . 10[) yields the first term missing in (|3.7p 

A 5 = -4m 4 7r 4 6. 

3.2. Construction of an approximate solution from a formal solution. For 

any order of approximation MeN, we now use the formal solution obtained above 
to construct a candidate for an approximate eigenvalue and eigenfunction on fijy. 
See Remark 1141 for a motivation of our matching procedure. In this subsection we 
still fix m and omit it from the notation. 

Choose cut-off functions Xf, Xg on as follows: Choose a smooth function x on 
R which equals one on (—00, 1/2) and zero on (3/4, 00). Then set Xfi x i V) = x( x /N) 
and Xg( x ,y) = 1 ~ XO)- 
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Set 

(3.11) / (M) = £ N-% 



M-l 



i=0 
M-l 

(3.12) g {M) (^y)=Y, N ' 1 9 J (^,y), 

3=0 
M-l 

(3.13) w {M) {x,y) = Y iV _ Vt(»)sin7r2/. 

W ( M ) describes the essential large x behavior of and the small £ behavior of 
g (M \ Set 

(3.14) f/ (M) =X// (M) +X 9 3 (M) -X/X 9 ^). 

That is, f/W is given by for x < 1/2, by g( M > for x > (3/4) JV, and by a 

smooth, appropriately scaled transition in between. Finally, set 

M-l 

(3.15) = 7T 2 + N ~" X i' 
Proposition 13. Denote 

V (M) = (A + A (M) )f/ (M) . 
For any m, M, p G N i/iere are constants c, C > smc/i f/iai /or aH AT we have 

(3.16) |A P T/ (M) | < C (log N/N) M , uniformly m Q N , 

(3.17) u (M) jA py(M) — o on d Q N ^ 

and 

(3-18) \\U^\\ LHnN) >cN^. 

All constants C, as well as cr x , are bounded in terms of max 0. 

Proof. The idea is to split up (|3.14p in two ways: First, as the / term plus g — w, 
which gives estimates on the order of (x/N) M and then as the g term plus / — w 
which gives estimates on the order of N~ M + e~ x . One of these is always bounded 
as in (|3~TC|) . 

Denote (x) = 1 + |x|. First, using the equations (f) for i < M we get 
(A + AW)/W= £ N-W hfi . 

i,l<M-l 
i+l>M 

Applying (f) again, we see that Afi is a linear combination of the fj, j < i, and 
by induction over p we get that A P (A + \(- M )~)f(. M ) is a linear combination of 
ft,..., /m-i, with coefhcients bounded by N~ M . This implies 

(3.19) A*(A + A(")) W /W = O(^L) 

uniformly in f2jv, since each = ©((a;) 1 ), and since derivatives of x/ are O^" 1 ) 
and only occur where x > N/2, where each N~ l fi, and therefore f^ M \ and its 
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derivatives of any order are uniformly bounded. Similarly, equation (g) and ()3.5j) 
give 

(A + A (M) )5 (M )= £ N-Wxjm, 

j,l<M-l 
j+l>M+2 

and using (g) again and boundedness of the gj we get 

(3.20) A p (A + \ (M) )xg9 (M) = 0(N- M - 2 ) for x > 1. 

Next, expanding gj(£, y) = Y^iLq 1 d \gj{0, y)€ l /l\ + £ M Rj(£, y), with Rj smooth, we 
get 

g (M)_ XfW (M) 

M-l f)l (n \ M-l M-l 

j,Z=0 i=0 i<i j=0 

Writing i = J +i in the first sum, we see from (Mu) that all terms of order at most 

(3-21) A^ g ( ff W-x/^ M) ) = 0(^). 

Finally, we have 

(3.22) K>Xf(f (M) -X 9 ™ (M) ) = 0(e- x ) for x > 1 



immediately from (I3.4[) and 

Writing £/( M ) = Xff (M) + Xg{g (M) ~ XfW (M} ), we get from (13J91) and (jiT2Tll 

(3.23) A p (A + A (M) )C7 (Af) = 

and writing JjW = x s .9 (M) + Xf(f (M) - XgW {M) ) we get from (|3~^0j) and (j3~2"2"j) 

(3.24) A P (A + A (M) )C/ (M) = 0(iV- M - 2 + e~ x ), x > I. 

Using (f3~2"3")l for x < MlogN and (|3~2"4"|) otherwise we obtain (|3~TB)) . 

The boundary conditions (|3 . 1T[) are also clear from the arguments above (note 
that x/j Xg only depend oni, so one only gets ^-derivatives of gj in the terms 
where the cut-offs are differentiated). 

(|3.18p follows immediately from the estimate | sin rmrx/N \ > sin7r/8 > 0, x E 
(N(1-^),N(1- g^)), which implies > l/10for these x and y S (1/4,3/4), 

for large N. □ 

Remark 14. Let us clarify our procedure of obtaining asymptotic eigenfunctions, 
by relating it to a simpler, 'compact' problem. First recall how one may obtain a 
smooth function u(s,t) on R 2 with given Taylor expansions u ~ Fi( s ) I a ^ 
t = and u ~ s^Gj^t)/ j\ at s = 0, at least up to a certain order M: First, such 
a m exists iff the mixed derivatives of u at (0, 0) obtained from the two expansions 
agree, that is if 

(3.25) 4Fi(0) = 4Gj{0) fora!H,j. 
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Calling this common value Wij and setting, for a given order of approximation M, 

p{M) = Z M-i tiFiWili g( m) = g-v^w/ii, «,w = z^Jot^m 

one may set 

U (M) = p(M) + Q(M) _ W (M) _ 

From 

^ - W ( M) = £ | - X) &v s n = o( S M ) 

1=0 l ' \ j=0 J ' J 

one sees that u( M ) - G( M > = 0(s M ) uniformly for (s,i) 

near zero, and similarly 

U ( M ) _ p(. M ) = 0{t M ) uniformly near zero, which was our goal. These estimates 
continue to hold if one formally differentiates both sides any number of times. 

This may be used to construct asymptotic solutions of partial differential equa- 
tions: Let P be a partial differential operator of b-type, i.e. a polynomial in sd s , 
tdt with smooth coefficients. Suppose one can determine the Fi and Gj so that 
PF( M ) = 0(t M ), PG {M) = 0(s M ), uniformly near zero, (which amounts to solving 
a recursive set of ordinary differential equations for the Fi and the Gj) then 

P U W = PG {M) + p (F (M) _ W (M) } = 0(s A/ } 

and similarly Pu ( - M '> = 0(t M ), so 

(3.26) Pu (M) = 0(min{s M , t M }) = 0((st) M/2 ). 

This relates to our problem as follows: We want to describe the eigenfunction u 
uniformly in x and N, that is, as a function on 

D={(N,x,y): (x, y) e Qn} C M 3 . 

In the sequel we suppress the y-dependence for simplicity. Our ansatz postulates 
that u has nice expansions in terms of smooth functions of x and x/N . This may 
be expressed as follows: Introduce new variables 



in the subset {x > 1} of D. Allowing the value iV = oo (i.e. adding fioo) and then 
s = 0, we get a compactification D of D, given by adding a point at infinity for 
each value of t S [0, 1] and y € [0, 1]. What we prove is that u extends to a function 
on D which is smooth in s and t up to s = 0, t = 0. The expansion (13. 2p may be 
rewritten 

s fixed, where Fi(s) = s l fi(-)il, 



t fixed, where Gj(t) = t J gj(t)jl. 

The matching conditions (Mu) are precisely the conditions (|3.25|) . with I = i — j. 
Note that Fi is smooth at s — by (f3~4|) . 

The cutoff functions in (I3.14|) must be introduced since F^ M ^ does not satisfy the 
boundary conditions at t = 1 and G^ M ^ does not satisfy the boundary conditions 
at s — > oo, i.e. at the left end of Qn- In (s,t) coordinates, the cut-offs are simply 
functions of s resp. t, and this motivates their choice in (|3. 14jl . 



i 

and (|3.3p becomes 
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We have d x = §^d s + §^dt = —s 2 d s + stdt, so the Laplacian is of b-type. The 
estimate f|3 . 1 6[) is actually stronger than what should be expected from (|3.26p (which 
gives 0{N~ M / 2 ) only) since in our problem the structure of P yields Gj = 0(t J ) 
andFM-«)W=0(e- 1 /*). 

3.3. Closeness to actual solution, proof of Theorem [1} 

Theorem 15. Denote the Dirichlet eigenvalues of —A on £In by fix < f/,2 < • • • , 
and denote the approximate rath eigenvalue constructed above by \nf\ and the 
approximate eigenfunction by Um 1 ^. Fix to. For sufficiently large N the first to 
eigenvalues on Qn are simple, and for each M 

(3.27) \^-X^ M) \=0(N- M ) forj = l,...,m. 

Furthermore, for j < m there is an eigenfunction Uj for the eigenvalue satisfying, 
for any a > 0, 

(3.28) S up\D a ( Uj -U^ M) )\=0(N- M ), j = l,...,m. 
The implied constants only depend on M,j,a and max0. 

Proof. Let i>i,i>2> • • • be an orthonormal basis of eigenfunctions on Qn, correspond- 
ing to the eigenvalues fa, fa, For fixed j £ {1, . . . , to} write 

(3.29) Uf ] =^ am . 

i 



Then ai = (£/} >/) (scalar product in L 2 {Q N )). For Vf> = (A + \y')U)j' we 



then obtain (V- M \vi) — (A^ M ' — ^i)ai, using integration by parts and U.ffl = 0, 



and then by induction (A p vj M \vi) — ^(X^ — /i;)a/ for all p > using (|3.17[) . 

since (I3~TCD implies ||APV 7 (M) || L 2 < CN- M+1 we get from Parseval's formula and 
(I3~18l) 

(3.30) JIl a '| 2 - cAr 

(3.31) I a » - mf < CN- 2M + 2 . 
i 

From (13301) and (l3~BT)l . with p = 0, we get \xf I] - fa\ < CN- M+1 / 2 for some I. 

Taking M > 4, we get that there is an eigenvalue flu in a CW~ 3 -neighborhood 
of it 2 + j 2 ir 2 /N 2 , for each j = 1, . . . , to. Since these neighborhoods are disjoint for 
N sufficiently large, we have < • • • < (J>i m , in particular pn m > /i m . On the 
other hand, comparing Q.^ to the larger domain [— C, iV] x [0, 1], we see that the 
TOth eigenvalue (counting multiplicity) of is at least tt 2 + m 2 ir 2 /(N + C) 2 — 
TT 2 +m 2 iT 2 /N 2 +0(N~ 3 ). Therefore,^™ > fa m . This implies that, for N sufficiently 
large, [ij — fii j for j = 1 . . . ,to , so the eigenvalues fa are all simple and \fa — 
A (M)| = ( N -M+i/2y R e pi ac i ng m by M + 1 and subtracting N~ M Xmj horn 

Xf 1+1) we obtain (l3~27l) . 

In particular, we have \fa — > cN~ 2 for Z 7^ j, and therefore we get from 

(13"3D) rf P \ai\ 2 < CN- 2M + 6 which means that 

\\A"(U^-a jVj )\\ LHQN) <CN- M + 3 . 
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By the Sobolev embedding theorem (applied to any unit width strip in Qn) we 
have ||r||o™(Ojv < C|| A p r|| L 2( njv) + C|MU 2 (n N ) for any function r on Q N , whenever 
2p > a + 1. Therefore, replacing M by M + 3 and then subtracting terms of order 
N~ l , i > M, on the left we get (|3 . 28[) with Uj — cijVj. □ 

Proof of Theorem^ (fTT2"j) follows with M = 5 from (f3~7]l . (|3~T5)l (where the index 
m was omitted in the notation) and (|3.27[) (for j = m). 

For the eigenfunction we first recall (|3.14[) , which gives for x > 1 (where \g ( x ) — 

1): 

(3-32) U iM )=g{ M )+Xf{f{ M ) _ w{ M )y 

From (|3"1^|) . (|3~5)) and we have 

(3.33) 5 (M) (x, y) = sin ™ir^- sin Try + 0(N~ 3 ), 

and from (13.22|) we have for x > 3 log iV 

(3.34) X/(/ (M) - w (M) ) = 0(iV~ 3 ). 

Clearly, the estimates Q3.33P and (|3.34[) may be differentiated any number of times. 
Therefore, (|3~28|) . (f3T32|) . ((3733]) and (|3~34|) give the eigenfunction estimate (fT73|) . 

Finally, from [0, iV] x [0,1] C Ojv C [— max <j>, AT] x [0,1] one has by domain 
monotonicity 

2 2 2 2 

2 TO Z 7T Z . ,„ N . , m Z 7T Z 



(iV + max0) 2 " ' TV 2 ' 

and combining this with (| 1 . 2[) one obtains < a{<f) < maxf □ 

4. Perturbation of the domain 
In this section we prove Theorem [2] 

First, we derive an alternative formula for a(<f>). (|2.2|) can be rewritten a(<f>) = 
2 J 9f2 U d(xsmTry)/dn, where d/dn denotes differentiation in direction of the out- 
ward unit normal n. Therefore, by applying Green's formula on Qq and using 
(A + n 2 )U — (A + 7r 2 ) (x sin Try) = we obtain 



a{4>) 

(4.1) 



„ f dU ■ 

2 / — — xsmTry dsiy) 

JdQo on 

(d x +4 > '(y)^v) u \{-4{y),y)4>{y) s'mwydy 



since n = (-1, -<j>')/y/l + ((j)') 2 and ds = y^TTW^dy. 

Now fix </> and denote by £l e N the domain Qn defined using e<f>, and by U e the 
associated function U from Proposition 

Note that Theorem [2] would follow from (14. 1|) (with replaced by e</>) if t/ = £/ £ 
could be replaced by U°(x,y) = x sin Try. Therefore, writing v e — U e — U we only 
need to show that 



(4.2) 



dn 



= O(e), 

i 2 (B e ) 



where i? e = {(— e<j){y), y) : y € [0, 1]} is the left boundary. Since v e BE = — e^sinTry, 
this follows from the following lemma. 
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Lemma 16. Suppose h is a function on dfloc, supported in B = {(—<fi(y),y) ■ y € 
[0, 1]}, and v solves 



then 



(A 



dv 



7T )v = on Qgo, 
v = h at dQco, 
v is bounded 



dn 



< C\\h\\ H i( B ) 



where C is bounded in terms of the Lipschitz constant of <fi. 

Proof. Write v = u + H , where H is an extension of h to Sloo, supported in fli, 
satisfying \\H \\ _f/i(Qi) < C||/i||#i(B). Then u satisfies the assumptions of Lemma[7] 
with w = -~(A + n 2 )H, so (pTT)) gives \\u\\ H i < \\w\\ H -i < \\H\\ H i and therefore 

(4-3) Mm^) <C\\h\\ m(B) . 

Next, we choose a smooth cut-off function x{ x )i equal to one in i < 1/2 and 
to zero in x > 3/4, and set v — %u. This satisfies Av = w, where w := —ir 2 v + 
2VxVv + (Ax)v, and V\qq 1 = h, so standard estimates give 



dv 



di 



LHB) 



< C(lklU 2 (ni) + WHh^b)) < C\\h\\ H i (B) 



using ||w||l2 < C\\v and (|4.3[) . Since v — v near B, this proves the lemma. □ 



5. Maximum set and nodal line 
Here we prove Theorem @] First, we obtain the following corollary of Theorem 

HI 

Corollary 17. Consider the eigen] 'unctions u\, ui on 
(a) // Mi assumes its maximum at a point (x, y) then 

N - a(4>) 



(5.1) 



■i = o(n 



(b) // {x,y) is an interior point offl^ with U2(x,y) = then 



(5.2) 



N - a(4>) 



\=0{N- 2 ). 



Proof. For shortness, we write a = a((f>). (a) First, from (|1.3p with a — and 
from (|1.4p it follows that, at a maximum (x, y), we must have x S [N/3, 2N/3] and 
sin7rj/ > 1/2, for large N. Next, we use that d x u\ = at a maximum. From (|1.3[) 
with a = (1,0), i.e. taking x-derivatives, we obtain after multiplication by N and 
division by sin7ry 

x + a ^, , T _ 2 \ 



N ■ 



0(i\T 



With x = ^-j^ + e the expression on the left becomes cos( : | + = — sin jf^, 

so from | sini| > \t\/2 for small t we get |e| = 0(N^ 1 ). 

(b) First, by integrating (|1.3|) with a = (0, 1) along the line from (x,0) to (x,y) 



we get \u m (x,y) 



. x+a(4>) 
N+a(4>) 



sin 7T?/ 1 < CyN 3 ; by a similar estimate in terms 
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Figure 2. Two domains with same height function but distinct 
location of the nodal line. 



of distance to the upper boundary y = 1 we obtain, after dividing by sin7ry, the 
improvement of |L 



(5.3) sup 

a;>31og N 



u m (x,y) . x + a(<t>) 

— : smmi- — - 

sm7ry Jy + a{<p) 



0(N- 3 ). 



If U2(x,y) — then, by Theorem 1 in (8], we have x E [JV/3, 2JV/3]. Therefore, we 
obtain from (|5.3p (with m = 2) 

sin27 r^±^ = 0(^ 3 ). 

N + a y ' 

As above this implies (|5.2|) . □ 

We consider domains Qn of the form (jl.ip which are convex, i.e. with a concave 
function cf>. By the corollary, to prove Theorem [4] it is enough to establish two 
concave functions <f>, <f> so that a{4>) ^ a(4>) and the corresponding domains SI at, fijv 
have the same projection and height function. 

Let (f>o(y) = 1/2 — \y — 1/2| and </>o(y) = 2//2- See Figure [2] Since (f> is the 
symmetric decreasing rearrangement of 4>q around the point 1/2, and since (sin7ry) 2 
is symmetric decreasing itself, we have L (j)o(y)(sm ny) 2 dy > f Q 4>o(y)(smny) 2 dy, 
so Theorem [5] implies a(e(f>o) ^ a(e^o) for some sufficiently small e > 0. Since 
the domains associated with 4> = e^o and 4> — e^o clearly have the same height 
function, the theorem is proved. 
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